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Abstract 

Acute  responses  to  intense  stressors  can  give  rise  to  post-traumatic  stress  disorder 
(PTSD).  PTSD  diagnostic  criteria  include  trauma  exposure  history  and  seif-reported  symp¬ 
toms.  Individuais  who  meet  PTSD  diagnostic  criteria  often  meet  criteria  for  additionai  psy¬ 
chiatric  diagnoses.  Biomarkers  promise  to  contribute  to  reiiable  phenotypes  of  PTSD  and 
comorbidities  by  linking  bioiogicai  system  alterations  to  behaviorai  symptoms.  Here  we 
have  anaiyzed  unbiased  piasma  metabolomics  and  other  stress  effects  in  a  mouse  modei 
with  behaviorai  features  of  PTSD.  In  this  model,  C57BL/6  mice  are  repeatediy  exposed  to  a 
trained  aggressor  mouse  (albino  SJL)  using  a  modified,  resident-intruder,  social  defeat  par¬ 
adigm.  Our  recent  studies  using  this  model  found  that  aggressor-exposed  mice  exhibited 
acute  stress  effects  including  changed  behaviors,  body  weight  gain,  increased  body  tem¬ 
perature,  as  weil  as  infiammatory  and  fibrotic  histopathoiogies  and  transcriptomic  changes 
of  heart  tissue.  Some  of  these  acute  stress  effects  persisted,  reminiscent  of  PTSD.  Here  we 
report  elevated  proteins  in  plasma  that  function  in  infiammation  and  responses  to  oxidative 
stress  and  damaged  tissue  at  24  hrs  post-stressor.  Additionai iy  at  this  acute  time  point, 
transcriptomic  analysis  indicated  liver  inflammation.  The  unbiased  metabolomics  analysis 
showed  altered  metabolites  in  plasma  at  24  hrs  that  only  partially  normalized  toward  control 
levels  after  stress-withdrawal  for  1 .5  or  4  wks.  In  particular,  gut-derived  metabolites  were  al¬ 
tered  at  24  hrs  post-stressor  and  remained  altered  up  to  4  wks  after  stress-withdrawal.  Also 
at  the  4  wk  time  point,  hyperlipidemia  and  suppressed  metabolites  of  amino  acids  and  car¬ 
bohydrates  in  plasma  coincided  with  transcriptomic  indicators  of  altered  liver  metabolism 
(activated  xenobiotic  and  lipid  metabolism).  Collectively,  these  system-wide  sequelae  to  re¬ 
peated  intense  stress  suggest  that  the  simultaneous  perturbed  functioning  of  multiple  organ 
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Introduction 

Post-traumatic  stress  disorder  (PTSD)  is  a  commonly  occurring,  debilitating  anxiety  disorder 
that  often  persists  for  many  years  and  is  frequently  associated  with  exposure  to  multiple  trau¬ 
mas  [1].  Significant  controversy  surrounds  the  diagnosis  of  PTSD  [2]  that  is  based  on  history 
of  exposure  to  extreme  traumatic  event(s)  and  reported  symptoms  including  persistent  re¬ 
experiencing  of  the  event(s),  persistent  avoidance  and  numbing  to  trauma-related  stimuli,  in¬ 
creased  arousal,  and  clinically  significant  distress  that  impairs  social,  occupational,  or  other  im¬ 
portant  areas  of  functioning  for  more  than  one  month  [3].  The  PTSD  lifetime  prevalence  rate 
from  the  National  Comorbidity  Survey  was  estimated  at  7.8%  among  the  general  population  of 
American  men  and  women  [4].  PTSD  prevalence  rates  are  higher  in  individuals  exposed  to  se¬ 
vere  traumas,  and  in  post-conflict  settings  ranged  from  16  to  37%  [5,6]  with  22%  of  289,328 
veterans  of  the  Iraq  and  Afghanistan  conflicts  entering  VA  health  care  from  2002  to  2008  after 
having  received  the  diagnosis  of  PTSD  [7]. 

Neurological/physiological  alterations  associated  with  PTSD  comprise  both  central  and  au¬ 
tonomic  nervous  systems,  and  include  hyperarousal  of  the  sympathetic  nervous  system  with 
sensitized  and  augmented  acoustic-startle  eye  hlink  reflex  and  sleep  ahnormalities,  as  well  as 
changes  in  noradrenergic,  hypothalamic-pituitary-adrenocortical,  serotonergic,  glutamatergic, 
thyroid,  endogenous  opioid,  and  other  systems  [6].  Brain  imaging  has  revealed  reduced  volume 
of  the  hippocampus  and  anterior  cingulate  cortex,  excessive  amygdala  activity,  and  reduced  ac¬ 
tivation  of  the  prefrontal  cortex. 

Psychometric  and  psychophysiologic  techniques  have  been  developed  to  assess  PTSD  [8]. 
Self-reporting  comprises  a  large  part  of  the  diagnosis  of  PTSD,  hut  it  is  prone  to  controversial 
diagnoses.  Additional  controversy  in  the  diagnosis  of  PTSD  relates  to  the  high  rate  of  comorbid 
conditions;  i.e.,  individuals  who  meet  the  PTSD  diagnostic  criteria  are  likely  to  meet  criteria  for 
one  or  more  additional  diagnoses  [6,9] .  As  such,  PTSD  may  he  categorized  into  multiple  di¬ 
mensions  with  traits  such  as  cognition,  mood,  and  social  interactions  existing  on  a  continuum 
from  normal  to  extreme  [10].  Ongoing  research  aims  to  continue  to  link  these  behavioral  di¬ 
mensions  with  underlying  molecular  pathologies  in  order  to  construct  valid,  measurable  traits 
for  PTSD  as  well  as  to  clarify  the  boundaries  and  overlap  with  other  mental  disorders,  toward 
improving  PTSD  prevention,  diagnosis  and  treatment  [10]. 

Animal  models  of  anxiety  disorders  vary  widely.  In  the  current  study,  we  used  a  previously 
validated  C57BL/6J  mouse  model  with  PTSD-like  features  [11].  Mice  are  stressed  by  exposures 
to  trained  aggressor  mice  for  5  or  10  6-hr  sessions  daily  in  a  cage-in-cage  environment  wherein 
the  subject  mouse  is  removed  from  the  smaller  cage  for  up  to  three  random  times  during  each 
6  hr  session  to  be  attacked  by  the  aggressor  for  one  minute  or  ten  strikes,  whichever  comes 
first.  In  previous  studies  with  this  model,  subject  mice  experienced  acute  stress  effects,  includ¬ 
ing  body  weight  gain,  increased  body  temperature,  cardiac  inflammation  and  transcriptomic 
changes,  and  behaviors  indicative  of  fear  and  anxiety.  Some  of  these  acute  stress  effects  per¬ 
sisted  following  1.5  to  6  weeks  of  stress-withdrawal,  reminiscent  of  PTSD.  Here  we  have  evalu¬ 
ated  additional  acute  and  chronic  stress  effects  in  this  model,  including  changes  in  the 
transcriptome  of  liver  and  proteins  and  metabolites  in  plasma.  In  particular,  transcriptomics 
indicated  liver  inflammation  at  24  hrs  post-stressor.  Unbiased  metabolomics  identified  altered 
metabolites  derived  from  the  gut  microbiota  as  early  as  24  hrs  post-stressor  that  remained  al¬ 
tered  up  to  4  wks.  Also,  at  this  chronic  time  point,  aggressor-exposed  mice  displayed 
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Figure  1 .  Study  design.  Mice  were  exposed  to  a  trained  aggressor  mouse  for  6  hrs  daily  for  5  or  1 0  days 
(black  rectangles).  Control  mice  were  exposed  on  the  same  schedule  and  under  the  same  conditions  but  with 
no  aggressor  mouse  present.  Terminal  bleeds  and  tissue  samples  (red  rectangles)  were  obtained  24  h  after 
the  5-  and  1 0-day  regimens  as  well  as  1 .5  weeks  after  the  5-day  and  4  weeks  after  the  1 0-day  stressor 
regimens.  Mice  were  singly  housed  throughout  the  experiment. 

doi:10.1371/journal.pone.0117092.g001 

metabolomic  markers  of  hyperlipidemia  and  transcriptomic  indicators  of  altered  xenobiotic 
and  lipid  metabolism  in  liver.  These  observations  contribute  to  defining  acute  stress  effects  and 
their  chronic  sequelae  following  repeated  exposures  to  intense  stressors. 

Results  and  Discussion 

Study  Design 

C57BL/6J  mice  were  analyzed  for  unbiased  plasma  metabolomics,  serum  proteins,  histopathol- 
ogy,  and  liver  transcriptomics  following  control  (Ctrl)  or  aggressor  exposure  (AggE)  for  6  hrs 
daily  for  either  5  or  10  days.  Twenty-four  hours  after  these  stressor  regimens,  as  well  as  1.5 
weeks  after  the  5-day  regimen,  and  4  weeks  after  the  10-day  regimen,  histopathology  was  per¬ 
formed  and  terminal  bleeds  and  tissues  were  obtained  for  analysis.  One  hour  prior  to  sample 
cohection,  during  a  partition  test  [11],  mice  were  exposed  to  the  potent  situational  cue  associat¬ 
ed  with  the  aggressor  exposures.  The  experimental  design  is  shown  in  Fig.  1. 

Plasma  metabolite  alterations  of  AggE  mice  at  24  hrs  post-stressor 
partially  normalized  after  stress- withdrawal  for  1 .5  or  4  wks 

The  unbiased  plasma  metabolomics  detected  330  named  metabolites  and  166  unidentified 
compounds  (SI  File).  For  the  named  metabolites,  the  pairwise  comparisons  of  their  levels  be¬ 
tween  Ctrl  and  AggE  were  associated  with  high  q-values,  likely  due  to  the  small  sample  size  of 
5  or  6  mice  per  group  and  the  known  highly  variable  response  of  even  inbred  mice  to  stress. 
However,  analysis  of  unadjusted  p-values  showed  trends  in  the  data.  Clearer  trends  were  iden¬ 
tified  by  analyzing  superpathways  and  subpathways  of  metabolites  constructed  based  on 
KEGG,  HMDB  and  the  internal  knowledge  base  of  Metabolon,  Inc.  (Durham,  NC).  Within 
many  of  these  pathways,  numbers  of  elevated  or  suppressed  metabolite  levels  in  AggE  vs.  Ctrl 
differed  from  the  null  hypothesis  of  equal  numbers  (binomial  test). 

In  pairwise  analysis,  AggE-altered  metabolites  were  more  numerous  at  24  hrs  after  the  last 
stress  session  as  compared  to  after  1.5  or  4  wks  of  stress-withdrawal  (Welch’s  t-test  and  2-way 
ANOVA):  40  and  14  metabolites  were  altered  at  24  hrs  and  1.5  wks,  respectively  (p  <  0.05), 
after  the  5-day  stress  regimen;  and,  37  and  20  metabolites  were  altered  at  24  hrs  and  4  wks,  re¬ 
spectively  {p<  0.05),  after  the  10-day  stress  regimen.  This  trend  suggests  that  the  acute  meta¬ 
bolic  perturbation  of  AggE  mice  partially  normalized  toward  the  Ctrl  level  after  the  mice  had 
been  withdrawn  from  the  stressor  for  1.5  and  4  wks. 
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Evidence  for  habituation  to  the  1 0-day  vs.  the  5-day  stressor  regimen 

Time- dependent  changes  in  plasma  metabolites  were  also  evident  in  both  Ctrl  and  AggE 
mice.  In  AggE  mice,  103  (5-day  regimen)  and  30  (10-day  regimen)  metabolite  levels  were 
significantly  different  between  the  acute  (24  hrs)  and  chronic  time  points  (1.5  or  4  wks). 
Similarly,  in  control  mice,  90  (5-day  regimen)  and  76  metabolites  differed  (10-day  regimen) 
between  the  acute  and  chronic  time  points  (SI  Table).  Thus,  greater  numbers  of  metabolites 
were  altered  with  time  after  the  5-day  regimen  as  compared  to  the  10-day  regimen,  despite  the 
10-day- regimen  mice  having  received  five  additional  AggE  sessions  and  having  been  withdrawn 
from  the  stress  for  2.5  wks  longer  than  the  5-day-regimen  mice.  This  may  indicate  that  the  mice 
receiving  the  10-day  regimen  habituated  to  the  stressor  between  the  5*  and  10*  daily  exposure 
(i.e.,  metabohtes  had  already  normalized  24  hrs  after  the  10*  stressor  session).  Consistent  with 
this  interpretation,  the  treatment  effect  (Ctrl  vs.  AggE,  acute  and  chronic  time  points  combined), 
the  time  effect  (acute  vs.  chronic;  Ctrl  and  AggE  combined),  and  the  treatment  and  time  interac¬ 
tion  (two-way  ANOVA),  were  aU  greater  after  the  5-day  compared  to  the  10-day  stress  regimen 
(SI  Table).  Thus,  there  were  fewer  metabohc  changes  24  hrs  after  the  10-day  stress  regimen  as 
compared  to  the  5 -day  stress  regimen,  suggesting  habituation  to  the  metabolic  effect  of 
the  stressor. 

Consistent  with  habituation,  the  percent  of  mice  that  jumped  in  response  to  the  aggressor 
or  fought  back  against  the  aggressor  dropped  between  the  4*  and  10*  days  of  the  stress 
regimen  (unpublished  data).  Additionally,  in  the  heart,  an  order  of  magnitude  fewer  AggE- 
vs.-Ctrl  transcripts  were  evident  at  24  hrs  after  the  10-day  than  after  the  5-day  regimen,  and 
inflammation  transcripts  peaked  on  the  third  day  of  the  AggE  regimen  and  then  fell  as  tissue 
remodeling  transcripts  rose,  suggesting  healing  despite  the  6  hr  daily  AggE  sessions  on  days  4 
to  10  [12]. 

As  a  further  indication  of  habituation,  in  the  liver  there  were  ~3-fold  more  AggE-vs.-Ctrl 
differentially  expressed  genes  (DEGs,  >  2-fold)  at  24  hrs  after  the  5-day  stressor  than  at  24  hrs 
after  the  10-day  stressor  (2169  vs.  668  for  5-  and  10-day  mice,  respectively;  Table  1).  Gene  set 
enrichment  analysis  of  these  DEGs  showed  the  downstream  functions  ‘biosynthesis  of  hydro¬ 
gen  peroxide’,  ‘accumulation  of  leukocytes’,  and  ‘activation  of  phagocytes’  as  the  top  three  and 
only  activated  functions  at  24  hrs  post-5-day  stressor  (category  ‘Physiological  System  Develop¬ 
ment  and  Function’,  (S2  Table).  Conversely,  at  24  hrs  after  the  10-day  stressor,  the  top  activat¬ 
ed  function  was  ‘binding  of  fibroblast  cell  lines’,  and  the  top  two  inhibited  functions  were 
‘quantity  of  neutrophils’  and  ‘quantity  of  phagocytes’.  These  results  are  consistent  with  those 
from  heart  tissue  (below  and  [12])  and  suggest  that  inflammation  of  livers  of  AggE  mice  at  24 
hrs  post-5-day  stressor  may  have  transitioned  to  fibroblast  infiltration  at  24  hrs  post- 10- 
day  stressor. 


Table  1.  Summary  of  differentially  expressed  genes  between  livers  of  Control  and  AggE  mice. 


5-Day  Stress  Regimen 

10-Day  Stress  Regimen 

24  hr 

1.5  wk 

24  hr 

4  wk 

Up-regulated  genes* 

842 

530 

383 

193 

Down-regulated  genes* 

1327 

522 

285 

139 

Totals 

2169 

1052 

668 

332 

Based  on  >  2  fold  change  of  probe  intensity  with  p-value  <  0.05. 
doi:10.1371/journal.pone.0117092.t001 
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Heart  histopathology  was  unique  to  AggE  mice;  lymphoplasmacytic 
infiltrates  of  other  organs  were  similar  in  Ctrl  and  AggE  mice  but  varied 
with  stressor  duration  and  post-stress  time 

Similar  to  the  heart  transcriptomics,  histopathology  of  heart  tissue  of  AggE  mice  at  24  hrs  after 
the  5-day  stress  regimen  appeared  to  be  a  higher  grade  than  after  the  10-day  regimen  [11]. 
Heart  inflammatory  histopathologies  of  the  5-day- AggE  mice  at  24  hrs  had  disappeared  by 
1.5  wks  when  fibrosis  or  fibroplasia  and  myocardial  degeneration  were  present.  In  10-day 
AggE  mice  at  the  24  hr  point,  much  more  fibrosis,  fibroplasia  and  degeneration  were  present 
compared  to  at  24  hrs  after  the  5-day  AggE,  indicating  the  injury  had  started  to  heal  by  the  1 1* 
day  after  the  start  of  the  10-day  regimen.  A  larger  percentage  of  5 -day  AggE  mice  at  1.5  wks 
than  the  10-day  AggE  mice  at  4  wks  showed  fibrosis  or  fibroplasia  and  myocardial  degenera¬ 
tion,  which  are  likely  remnants  of  the  injury  that  persisted  at  1.5  wks  but  was  evidently  healed 
by  4  wks  (SI  Fig.). 

Acute,  intense  psychological  stressors  such  as  inter-male  aggressive  encounters  or  catechol¬ 
amine  administration  have  been  shown  to  induce  organ  tissue  injury  [13,14].  Upon  fighting. 
Epidermal  Growth  Factor  (EGF)  released  from  submandibular  salivary  glands  protected 
against  heart  injury,  but  not  against  liver  injury  [15].  However,  twice  daily  adrenergic  stimula¬ 
tion  decreased  the  EGF  content  of  submandibular  glands  by  almost  90%,  and  further  adrener¬ 
gic  stimulation  did  not  increase  plasma  EGF  levels  [15].  This  result  suggests  protective 
mechanisms  against  acute  stress  can  become  exhausted  in  response  to  chronic  stress.  In  our 
model,  the  repeated,  intense,  5-  and  10-day  AggE  sessions  likely  exhausted  mechanisms  for 
protecting  heart  tissue  from  intense  adrenergic  stimulation. 

In  contrast  to  the  cardiac  lesions  found  only  in  AggE  mice,  lymphoplasmacytic  infiltrates  of 
liver,  gallbladder,  kidney,  and  stomach  were  of  similar  grades  in  Ctrl  and  AggE  mice.  This  sug¬ 
gests  the  possibility  that  the  Ctrl  condition,  which  reproduced  the  AggE  condition  with  the  ex¬ 
ception  that  the  aggressor  was  absent,  was  a  significant  stressor,  likely  involving  social  isolation 
stress.  Although  similar  in  Ctrl  and  AggE  mice,  immune  cell  infiltrates  of  liver,  gallbladder, 
kidney,  and  stomach  differed  from  each  other  in  trends  of  rate  of  development  and  the  impact 
of  the  stressor  duration  (S2  Fig.).  Kidney  lymphoplasmacytic  interstitial  infiltrates  apparently 
worsened  with  time,  while  liver  lymphoplasmacytic  infiltrates  lessened  during  stress-withdraw¬ 
al.  For  organs  with  infiltrates  that  increased  during  stress-withdrawal,  the  trends  suggest  the 
possibility  that  the  10-day  regimen  might  have  been  more  damaging  than  the  5-day  regimen, 
and/or  the  alterations  could  have  been  cumulative  or  developed  more  fuUy  by  4  wks  than  by 
1.5  wks. 

Elevated  plasma  proteins  at  24  hrs  post-stressor  indicate  damaged 
tissue  in  AggE  mice 

In  plasma,  we  detected  41  proteins  in  the  Rodent  Multi- Analyte  Profile  v.2.0  Antigens  (Myri- 
ad-RBM)  panel  at  24  hrs  after  the  5-  and  10-day  stressors.  This  analysis  showed  that  haptoglo¬ 
bin  and  myeloperoxidase  and  serum  amyloid  P-component  (False  Discovery  Rate  p<0.05 
[16])  were  significantly  elevated  in  plasma  of  AggE  over  Ctrl  mice  at  24  hrs  for  the  10-day- 
exposed  mice  as  well  as  for  the  combined  group  of  5-day-  and  10-day-exposed  mice  (S2  File). 
These  proteins  are  involved  in  responses  to  oxidative  stress,  damaged  tissue,  and  inflammation. 
Myeloperoxidase  is  stored  in  azurophilic  granules  of  neutrophils,  oxidizes  chloride  ions  to  the 
potent  bactericidal  oxidant  hypochlorous  acid  [17],  and  has  been  widely  used  as  marker  of  in¬ 
flammation  in  acute  and  chronic  conditions,  including  cardiovascular  disease  [18].  Haptoglo¬ 
bin  is  an  acute  phase  protein,  increased  during  inflammation  [19],  which  has  been  observed  to 
rise  in  response  to  a  number  of  psychological  stresses  as  well  as  depression  [20,21].  Serum 
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amyloid  P-component  (SAP)  is  a  member  of  the  highly  conserved  pentraxin  superfamily  along 
with  the  classical  acute-phase  reactant,  C-reactive  protein  (CRP);  both  SAP  and  CRP  transcrip¬ 
tion  are  activated  in  response  to  endoplasmic  reticulum  stress  [22].  SAP  has  anti-fibrotic  effects 
[23]  and  binds  molecular  arrays  such  as  DNA,  chromatin,  histones,  and  phosphoethanolamine- 
containing  membranes,  suggesting  it  functions  in  the  clearance  of  late  apoptotic  cells  [24]. 

Additionally,  tissue  inhibitor  of  metalloproteinases  1  (TIMP-1)  was  elevated  in  10-day 
AggE  over  Ctrl  mice  (FDR  corrected,  p  =  0.038),  and  matrix  metalloproteinase-9  (MMP-9)  was 
elevated  in  the  5-day-exposed  over  the  10-day-exposed  mice  (combined  Ctrl  and  AggE 
groups).  MMP-9  and  matrix  metalloproteinase  family  members  are  involved  in  the  degrada¬ 
tion  of  extracellular  matrix  in  processes  such  as  wound  healing,  and  TlMP-1  is  natural  inhibi¬ 
tor  of  MMP-9,  and  the  two  proteins  are  markers  consistently  indicted  in  cardiovascular  disease 
development  and  prognosis  [25,26].  MMP-9  levels  have  been  associated  with  direct  effects  of 
cortisol  and  norepinephrine  (references  in  [27]),  and  independently  predicted  new  cardiac 
events  in  patients  with  stable  coronary  disease  [28],  as  well  as  having  been  independently  asso¬ 
ciated  with  psychosocial  factors  in  a  normal  middle-aged  population  [27].  Additionally,  in  our 
model,  in  heart  tissues,  MMP-2  and  TIMP-1  transcripts  were  elevated  in  AggE  mice  at  24  hrs 
after  the  5-day  and  10-day  stress  regimens. 

In  contrast,  granulocyte  chemotactic  protein-2  (GCP-2)  and  macrophage-derived  chemo- 
kine  (MDC)  were  elevated  in  Ctrl  relative  to  AggE  mice  at  24  hrs,  which  may  indicate  a  differ¬ 
ent  trajectory  of  inflammation  resulting  from  the  stress  of  the  Ctrl  condition. 


Gut  microbiota-derived  metabolites  in  plasma  of  AggE  mice  were  altered 
at  both  acute  (24  hrs)  and  chronic  (1 .5&4  wks)  time  points 

To  explore  metabolic  alterations  occurring  at  24  hrs  that  potentially  persisted  out  to  1.5  and  4 
wks,  we  used  random  forests  (RE)  to  classify  AggE  vs.  Ctrl  at  24  hrs  and  1.5&4  wks  (combined 
5-day  and  10-day  groups).  The  resulting  lists  of  the  top  10%  of  metabolites  with  the  greatest 
‘mean  decrease  accuracy’  (MDA)  at  the  acute  and  chronic  time  points  were  analyzed  for  over¬ 
laps  [29].  Nine  metabolites  overlapped  between  the  two  lists  (Table  2,p-  0.002,  Fisher’s  exact 
test).  Five  of  the  9  overlapping  metabolites  were  gut-derived  metabolites:  phenylpropionylgly- 
cine,  phenol  sulfate,  hippurate,  3-phenylpropionate,  and  p-cresol  sulfate.  Three  of  these — 3- 
phenylpropionate,  phenylpropionylglycine,  and  hippurate — were  elevated  in  AggE  mice  at 
both  acute  and  chronic  time  points  post- 10-day  regimen  (t-test,p  <  0.05),  and  the  acute  and 
chronic  time  points  of  the  combined  5-  and  10-day  group  (t-test,p  <  0.05).  These  three  metab¬ 
olites  are  among  the  most  elevated  metabolites  in  AggE  mice,  being  on  average  6.7-fold  elevat¬ 
ed  at  24  hrs  and  2.8-fold  elevated  at  1.5&4wks. 

The  metabolite  3-phenylpropionate  is  a  known  metabolic  product  of  anaerobic  bacteria 
[30].  Both  hippurate  and  phenylpropionylglycine  are  derived  from  3-phenylpropionate. 
Phenylpropionylglycine  is  the  glycine  conjugate  of  3-phenylpropionate.  Hippurate  is  a  glycine 
conjugate  of  benzoate,  derived  from  3-phenylpropionate  via  catalysis  by  Medium-chain 
acyl-CoA  dehydrogenase  (MCAD)  [31,32].  Hippurate  has  been  found  to  be  17-fold  lower  in 
plasma  of  germ-free  mice  than  conventional  mice,  and  phenylpropionylglycine  was  detected 
exclusively  in  conventional  samples  [30].  Phenol  sulfate  and  p-cresol  sulfate  have  been  shown 
to  be  reduced  several  hundred  to  1000-fold  in  germ-free  mice  [30]. 

The  remaining  four  metabolites  identified  by  RF  that  overlapped  between  24  hrs  and  1.5&4 
wks  were  apparently  diverse.  Notably,  2’-deoxycytidine  was  the  best  single  classifier  of  AggE 
vs.  Ctrl  mice  in  the  combined  5-  and  10-day  groups  at  24  hrs  using  either  distance-dependent 
K-nearest  neighbor  analysis  (sensitivity  =  90.9%,  specificity  =  81.8%)  or  rule-based  machine 
learning. 
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Table  2.  Metabolites  that  classify  AggE  vs.  Ctrl  mice  at  both  acute  and  chronic  time  points. 


24  hrs 

1.5&4  wks 

MDA 

FC 

p-value 

MDA 

FC 

p-value 

Phenylpropionylglycine* 

14.3 

10 

0.005 

3.5 

3.95 

0.024 

phenol  sulfate* 

14.2 

1.27 

0.090 

5.9 

0.77 

0.143 

Hippurate* 

12.6 

8.07 

0.008 

3 

2.52 

0.025 

3-phenylpropionate  (hydrocinnamate)* 

7.9 

2.16 

0.021 

7.3 

1.97 

0.015 

p-cresol  sulfate* 

7.1 

0.69 

0.244 

5.2 

0.77 

0.054 

2'-deoxycytidine 

22 

1.64 

0.001 

2.1 

1.13 

0.152 

palmitoyl  sphingomyelin  (d18:1/16:0) 

13.5 

1.26 

0.007 

15 

1.11 

0.015 

creatinine 

7.8 

0.76 

0.022 

8.2 

0.88 

0.159 

2-linoleoylglycerophosphocholine 

6.8 

0.83 

0.183 

6.1 

0.91 

0.277 

MDA — ‘mean  decrease  accuracy’  from  random  forests  analysis 
FC — fold  change 

*metabolites  derived  from  gut  microbiota 


doi:1 0.1 371/journal. pone.01 1 7092.t002 


Gut  microbiota-derived  and  other  metabolites  in  plasma  altered  at  either 
acute  or  chronic  time  points 

In  addition  to  the  gut-derived  metabolites  identified  by  RF,  pairwise  and  correlation  analyses 
revealed  additional  gut-derived  metabolites  that  were  elevated  in  AggE  at  24  hrs.  Two  of  six 
metabolites  that  were  “perfect  classifiers”  of  AggE  vs.  Ctrl  mice  (i.e.,  levels  were  all  higher  or 
lower  in  AggE  mice)  at  24  hrs  post-5-day  stressor  are  metabolites  of  the  gut  microbiota;  2-(4- 
hydroxyphenyl)propionate  and  indolelactate.  Both  were  elevated  in  AggE  vs.  Ctrl  mice  (com¬ 
bined  5-  and  10-day  groups).  2-4-hydroxyphenyl)propionate  is  a  fecal  bacteria  fermentation 
product  of  tyrosine  [33,34],  which  increased  in  feces  of  antibiotic-treated  mice  [35],  and  it  also 
positively  correlated  in  our  dataset  with  indolelactate  (r  =  0.86),  which  is  the  other  elevated  per¬ 
fect  classifier.  Indolelactate  is  a  tryptophan  metabolite  produced  by  Bifidobacterium  species 
[36],  and  was  found  to  be  at  least  10-fold  more  abundant  in  urine  of  germ-free  mice  colonized 
with  Bifidobacterium  longum  than  in  germ-free  mice  or  mice  colonized  with  Bacteroides  the- 
taiotaomicron  [37].  It  was  also  7-fold  reduced  in  serum  in  a  model  of  inflammatory  bowel  dis¬ 
ease  [38].  Additionally,  2-(4-hydroxyphenyl)propionate  correlated  with  indolepropionate  (r  = 
0.71).  Indolepropionate  has  been  identified  in  the  plasma  of  conventional  mice  but  not  in 
germ-free  mice  [30]. 

Of  the  top  10%,  or  33,  highest  MDA  metabolites,  24  did  not  overlap  between  the  24  hr  and 
1.5&4  wk  time  points.  Unique  to  the  24  hr  time  points  were  the  additional  gut-derived  metabo¬ 
lites  indolelactate  and  2 -(4-hydroxyphenyl) -propionate  (also  identified  as  perfecl  classifiers), 
and  one  of  their  correlating  gut-derived  metabolites,  indolepropionate.  Another  gut-derived 
metabolite  altered  at  24  hrs  was  phenyllactate  [39].  Thus,  these  gut-derived  metabolites  appear 
altered  at  the  acute  time,  in  addition  to  the  gut-derived  metabolites  altered  at  both  the  acute 
and  chronic  times.  Additionally,  two  fibrinogen  cleavage  peptides  were  highly  ranked  unique 
RE  classifiers  of  AggE  mice  at  the  24  hr  time  point.  Fibrinogen  cleavage  products  have  been  as¬ 
sociated  with  inflammation  [40]. 

The  RF- classifying  metabolites  unique  to  1.5&4wks  included  three  plant-derived  biochemi¬ 
cals,  campesterol  [41],  piperidine  [42]  and  daidzein  [43].  Additionally,  seven  lysophospholipids 
were  among  the  RF  top-scoring  metabolites  at  1.5&4  wks  and  they  were  all  elevated.  These  re¬ 
sults  suggest  that  AggE  stress  may  have  damaged  the  GI  tract  [44].  In  addition  to  a  likely 
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altered  gut  microbiota,  other  changes  could  possibly  include  altered  intestinal  absorption  and 
altered  liver  metabolism.  The  GI  tract  responds  to  stress  and  mediators  of  stress  (including  glu¬ 
cocorticoids  and  catecholamines)  with  disturbed  motor  activity,  ulcerations,  enteritis,  and  al¬ 
tered  transepithelial  ion  transport  [45-47],  all  of  which  can  affect  the  composition  of  the 
microbiota  [48].  Also,  gut  bacteria  can  respond  directly  to  catecholamines  with  changes  in 
growth  and  virulence  gene  expression  [45].  Alterations  in  intestinal  bacterial-host  interactions 
have  been  suggested  to  cause  liver  injury  [49]. 


Metabolites  of  the  major  metabolic  fuels  in  plasma  were  altered  at  acute 
and  chronic  time  points 

The  major  metabolic  fuels,  carbohydrates,  amino  acids,  and  lipids  differed  in  plasma  levels  de¬ 
pendent  upon  the  stressor  duration  and  the  length  of  the  stress-withdrawal  period  (Fig.  2).  At 
24  hrs  after  the  5-  and  10-day  stressors,  81%  and  76%  of  the  21  carbohydrate  metabolites  were 
higher  in  plasma  of  AggE  mice  (binomial  test,p  =  0.003,  post-5-day  stressor,  andp  =  0.013, 
post-lO-day  stressor).  But  after  withdrawal  from  the  stress  for  1.5  and  4  wks,  only  57%  and 
29%  of  plasma  carbohydrates  were  elevated,  fewer  than  at  24  hrs  (test  of  equal  proportions:  p  - 
0.182  andp  =  0.005,  respectively;  and p  -  0.002  for  the  combined  5-day  and  10-day  group). 

Acutely  elevated  carbohydrate  metabolites  in  plasma  of  AggE  mice  could  possibly  reflect  the 
action  of  the  primary  stress  mediators,  epinephrine  and  corticosterone,  on  fuel  storage  and 
availability.  Corticosterone  mobilizes  amino  acids  mainly  from  muscle,  making  more  available 
in  the  plasma  to  enter  gluconeogenesis  in  the  liver  to  promote  glucose  and  glycogen  synthesis 
[50].  Thus,  more  glycogen  is  made  available  for  glycogenolysis  stimulated  by  epinephrine, 
which  can  release  within  minutes  large  amounts  of  glucose  into  the  blood  [50].  The  rise  in 
blood  glucose  in  turn  stimulates  insulin  secretion,  but  high  levels  of  glucocorticoid  make  skele¬ 
tal  muscle  and  adipose  tissue  resistant  to  insulin’s  stimulatory  effects  on  glucose  uptake  and 
utilization  [50].  Consistent  with  this  scenario,  at  the  acute  time  point,  glucose  level  showed  a 
trend  to  be  elevated  1.2-fold  in  the  5-day  24  hr  group  (p  <  0.1)  and  1.18-fold  in  the  combined 
5-day  and  10-day  24  hr  group  {p  <  0.1).  Also  possibly  consistent  with  this  trend,  2-hydroxybu- 
tyrate,  an  indicator  of  insulin  resistance  and  impaired  glucose  regulation  [51-53],  was  1.74- 
fold  elevated  in  the  combined  group  of  5-  and  10-day  AggE  mice  at  24  hrs  {p  -  0.007).  Howev¬ 
er,  changes  in  2-hydroxybutyrate  have  also  been  related  to  the  oxidative  stress  response  in  the 
liver  and  increased  glutathione  synthesis. 

Conversely,  the  number  of  amino  acid  metabolites  elevated  in  AggE  mice  at  24  hrs  post-5- 
day  stressor  was  not  greater,  but  more  of  them  were  elevated  at  24  hrs  after  the  10-day  stressor 


90 


24  hrs  LSwAs  24  hrs  4  whs 

S- Day  Stressor  lO-Oay  Stressor 

■  Carbohydrates  ■  Ammo  Acids  ■  Lipids 


Figure  2.  Percent  of  metabolites  of  metabolic  fuels  higher  in  AggE  vs.  Ctrl  mice.  The  Carbohydrate, 
Amino  acid  and  Lipid  superpathways  are  comprised  of  21 , 98,  and  1 34  metabolites,  respectively.  Asterisks 
indicate  p<0.05  (*)  and  p<0.01  (**)  for  obtaining  the  percent  elevated  metabolites  shown  (exact  binomial 
test). 

doi:10.1371/journal.pone.0117092.g002 
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(60%  of  98  metabolites,/)  =  0.017).  Of  the  major  metabolic  fuels,  only  amino  acid  metabolites 
appeared  to  differ  between  the  5-day-24  hr  and  10-day- 24  hr  groups  (proportion  test,p  = 
0.063),  coincident  with  habituation.  More  amino  acid  metabolites  were  also  elevated  in  AggE 
mice  at  1.5  wks,  and  significantly  more  than  at  24  hrs  {p  -  0.001).  The  fewer  plasma  amino 
acid  metabolites  of  AggE  mice  at  24  hours  after  the  5-day  regimen  may  have  possibly  resulted 
from  higher  glucocorticoid  signaling,  resulting  in  depleted  muscle  capacity  for  mobilizing 
amino  acids. 

At  4  wks  post-stressor,  the  great  majority  of  amino  acid  metabolites,  as  well  as  carbohydrate 
metabolites,  were  lower  in  the  plasma  of  AggE  vs.  Ctrl  mice.  Only  29%  of  carbohydrate  and 
16%  of  amino  acid  metabolites  were  higher  in  AggE  than  in  Ctrl  mice  at  this  time  point,  while 
68%  of  lipid  metabolites  were  elevated  in  AggE. 

A  possible  explanation  for  these  trends  is  that  at  the  4  wk  time  point,  when  the  allostatic 
load  of  AggE  mice  was  evidently  lowest,  there  may  have  been  a  return  toward  baseline  of  the 
sensitivity  of  hypothalamic-pituitary  axis  (HP A),  as  well  as  its  effectors  muscle  and  liver.  In 
this  scenario,  the  elevated  lipid  metabolites  may  have  resulted  in  part  from  the  situational  cue 
of  the  aggressor  during  the  partition  test,  which  occurred  less  than  one  hour  prior  to  the  termi¬ 
nal  bleed.  This  likely  elicited  a  fear-related  memory  that  activated  the  sympathetic-adrenome- 
duUary  axis,  thus  producing  epinephrine,  which  activates  hormone-sensitive  lipase  to  mobilize 
lipids  from  adipose  tissue.  In  contrast,  glucocorticoid  genomic  programing  of  muscle  and  liver 
at  4  wks  had  likely  returned  nearer  to  baseline,  and  the  fear-related  memory  only  one  hour  be¬ 
fore  the  terminal  bleed  may  have  been  too  short  for  slow  glucocorticoid  genomic  changes  to 
have  taken  fuU  effect.  The  results  of  this  could  have  been  lower  amino  acid  mobilization  and 
gluconeogenesis,  which  would  explain  the  low  numbers  of  amino  acid  and  carbohydrate  me¬ 
tabolites  in  the  plasma  of  AggE  mice  at  4  wks  relative  to  the  other  time  points  when  allostatic 
load  was  evidently  higher. 

Consistent  with  this  scenario,  plasma  corticosterone  trended  higher  at  4  wks  relative  to  the 
earlier  time  points.  In  our  unbiased  metabolomic  analysis  (semi-quantitative  LC-MS/MS  anal¬ 
ysis),  corticosterone  trended  lower  in  AggE  than  in  Ctrl  mice  at  the  two  24  hr  points  and  at 
1.5  wks,  but  trended  higher  in  the  AggE  mice  at  4  wks.  Following  up  on  this  result,  we  per¬ 
formed  a  corticosterone  ELISA  on  plasma  samples  from  a  replicate  set  of  mice  groups  and  got 
the  same  result;  corticosterone  trended  lower  in  AggE  mice  at  the  two  24-hour  points  and  at 
1.5  weeks,  but  trended  higher  in  AggE  mice  at  4  wks.  Thus,  two  different  assays  performed  on 
different  plasma  samples  resulted  in  the  same  trend  of  relative  corticosterone  levels  between 
AggE  and  Ctrl  at  four  time  points  {p  -  0.06  for  this  binomial  distribution). 

One  possible  explanation  for  these  trends  is  that  at  both  24  hr  time  points  and  the  1.5  wks 
time  point,  the  HPA  had  already  become  desensitized,  potentially  involving  attenuated  release 
of  pituitary  ACTH  and/or  reduced  adrenal  responsiveness  [54].  However,  at  these  earlier  times 
of  higher  allostatic  load,  glucocorticoid-programing  of  muscle  and  liver  may  have  been  already 
put  them  into  a  heightened  state  for  mobilizing  amino  acids  and  producing  glycogen.  But  these 
heightened  states  could  have  diminished  by  4  wks,  although  HPA  responsiveness  may  have 
been  restored.  This  could  possibly  explain  the  elevated  corticosterone  one  hour  after  exposure 
to  the  situational  cue  of  the  AggE,  despite  carbohydrate  and  amino  acid  metabolites  being 
lower  in  AggE  mice  at  this  time. 

Another  model  based  on  the  social  defeat  paradigm  (in  NMRI  mice)  found  that  repeated 
stress  increased  core  body  temperature  (as  in  our  model)  and  after  one  social  defeat,  corticoste¬ 
rone  levels  peaked  at  30  min  [55].  By  one  day  after  the  single  social  defeat,  corticosterone  level 
was  also  increased,  but  not  after  the  4*  or  8*  defeat  over  2  wks,  while  it  was  dramatically  in¬ 
creased  24  hrs  after  the  12™  defeat  after  3  weeks  of  experimentation  [55].  Also  in  this  study, 
mice  after  the  12*  defeat  that  were  withdrawn  from  the  stress  for  9  days  were  sensitized  to  a 
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novel  stressor.  The  trend  we  observed  is  similar  to  this  and  suggests  the  possibility  that  the 
HPA  axis  was  de-sensitized  at  24  hrs  and  1.5  wks  and  may  have  become  re-sensitized  by  4  wks. 

Lipid  metabolites  in  plasma  of  AggE  mouse  were  elevated  at  all  the  time  points.  This  effect 
could  have  been  due  to  the  rapid  activation  of  hormone-sensitive  lipase  by  adrenalin  in 
adipose  tissue  to  increase  fat  utilization  for  energy  during  the  stress  of  the  fear  memory  and 
afterwards.  Elevated  lipid  metabolites  could  also  result  from  chronic  peripheral  effects,  in¬ 
cluding  damage  to  and  reprogramming  of  the  organs  that  absorb,  mobilize,  and  utilize  the 
major  fuels  [56]. 


Specific  metabolite  subpathways  altered  by  AggE  stress 

To  assess  possible  coordinate  regulation  of  metabolites  related  through  chemical  pathway  or 
structure,  we  used  a  Manhattan  distance  method  (see  Methods  section).  This  method  identified 
15  subpathways  (of  the  65  subpathways  in  our  data  set)  whose  metabolites  were  coordinately 
altered  across  the  8  comparisons  of  AggE  vs.  Ctrl  and  24  hrs  vs.  1.5  or  4  wks.  These  15  subpath¬ 
ways  were  hierarchically  clustered  based  on  their  metabolites’  distance  means.  This  produced 
two  large  clusters.  One  consists  of  9  subpathways  of  the  lipid  superpathway.  The  other  consists 
of  3  subpathways  under  the  Amino  Acid  superpathway  and  one  subpathway  under  each  of  the 
Peptide,  Carbohydrate  and  Xenobiotic  superpathways  (Fig.  3).  The  most  distinct  differences 
between  these  two  clusters  are  the  suppressed  metabolites  of  the  non-lipid  subpathways  and 
the  elevated  metabolites  of  the  lipid  subpathways  at  4  wks. 

We  also  evaluated  all  of  the  65  subpathways,  including  the  45  subpathways  without  signifi¬ 
cantly  co-regulated  metabolites,  to  identify  subpathways  with  more  AggE  vs.  Ctrl  metabolites 
elevated  or  suppressed  at  any  time  point  (exact  binomial  test).  The  most  marked  change  was  at 
4  wks  (Fig.  4)  when  many  of  the  altered  pathways  were  the  same  ones  found  with  the  Manhat¬ 
tan  distance  method.  Additionally,  the  exact  binomial  test  identified  the  subpathways  ‘Glycine, 
serine  and  threonine  metabolism’,  ‘Urea  cycle;  arginine-,  proline-,  metabolism’,  and  ‘Krebs 
cycle’,  as  being  suppressed  in  AggE  vs.  Ctrl  mice  at  4  wks  after  the  10-day  stressor.  The  sub¬ 
pathways  identified  by  both  methods  are  shown  overlaid  onto  the  network  of  their  structural 
and  pathway  relationships,  which  also  depicts  the  relative  levels  of  their  metabolites  in  plasma 
of  AggE  vs.  Ctrl  mice  (Fig.  5). 

Overall,  at  4  wks  after  the  stressor,  the  elevated  lipids  of  AggE  mice  plasma  included  several 
structural  classes  of  fatty  acids  (medium  and  long  chain,  essential,  and  monohydroxy),  lyso- 
phospholipids,  carnitine  metabolites  and  bile  acids.  Conversely,  suppressed  subpathways  in¬ 
cluded  ‘Glycolysis,  gluconeogenesis,  pyruvate  metabolism’,  ‘Krebs  cycle’  and  several 
subpathways  of  amino  acid  metabolism.  This  analysis  is  consistent  with  the  superpathway 
analysis  (major  fuel  metabolites)  discussed  above,  in  which  carbohydrates  and  amino  acids 
were  suppressed  and  lipids  were  elevated  in  AggE  vs.  Ctrl  mice  at  the  4  wk  time  point. 

The  relationships  of  these  altered  subpathways  to  each  other,  and  to  the  stressor  durations 
and  stress-withdrawal  periods  are  complex,  but  hypotheses  can  be  inferred.  For  example,  the 
five  conjugated  bile  acids  in  plasma  were  all  lower  {p  -  0.031,  binomial  test)  in  5-day  AggE 
mice  at  24  hrs  (by  1.32-fold  on  average),  which  could  reflect  higher  sympathetic  nervous  sys¬ 
tem  activity  thus  causing  changes  in  the  gut  and  liver  at  this  time  of  high  aUostatic  load.  In  con¬ 
trast,  after  a  1.5  wk  withdrawal  from  stress,  all  five  conjugated  bile  acids  were  higher  in  AggE 
mice  vs.  Ctrl  mice  (5.9-fold  average)  as  were  the  three  unconjugated  amino  acids  (i.e.,  8  of  8 
were  elevated  p  =  0.004).  Also,  in  the  combined  group  of  5-day  and  10-day  mice,  total  bile 
acids  in  AggE  plasma  was  1.15-fold  lower  (6  of  8  lower)  at  24  hrs  and  3.12-fold  higher  (6  of  8 
higher)  than  in  Ctrl  plasma  at  1.5&4  wks.  Elevated  total  serum  bile  acids  is  an  indicator  of  liver 
disease  and  can  reflect  reduced  hepatic  clearance  of  bile  acids  from  the  hepatic  portal  blood. 
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Figure  3.  Dendrogram  of  15  subpathways  whose  metabolites’  plasma  levels  changed  coordinately. 

Shown  is  the  hierarchical  clustering  of  subpathways  by  mean  Manhattan  distance  score,  which  describes  the 
similarity  in  the  direction  of  change  (higher  or  lower)  of  metabolites  of  a  subpathway  in  the  eight  comparisons 
of  AggE  vs.  Ctrl  and  24  hrs  vs.  1 .5  or  4  wks.  Also  shown  for  each  subpathway  is  the  percentage  of 
metabolites  higher  (red)  or  lower  (green)  in  AggE  vs.  Ctrl  mice  at  24  hrs  and  4  wks  after  the  10-day  stressor. 
The  distance  means  for  the  lipid  and  non-lipid  branches  of  the  dendrogram  were  3.1  (p  =  1 E-04)  and  3.07  (p  = 
1  E-04).  Asterisks  indicate  the  subpathways  also  determined  by  exact  binomial  test  to  have  more  metabolites 
lower  or  higher  in  AggE  vs.  Ctrl. 

cloi:10.1371/journal.pone.0117092.g003 


Additionally,  bile  acid  metabolism  can  be  regulated  by  the  gut  microbiota,  and  the  bile  acid 
tauro-muricholic  acid  is  an  antagonist  of  the  farnesoid  X  receptor  (FXR)  [57]. 

In  liver,  inflammation  at  24  hrs  and  metabolic  alterations  at  4  wks  are 
indicated  by  transcriptome  analysis 

Although  lymphoplasmacytic  infiltrates  of  livers  of  Ctrl  and  AggE  mice  showed  similar  levels 
(S2  Fig.),  sub-histological  alterations  could  have  been  present.  We  therefore  analyzed  differen¬ 
tially  expressed  genes  (DEGs)  of  liver  for  enrichment  in  the  ‘inflammatory  response’  function 
annotation.  More  AggE  vs.  Ctrl  DEGs  (up-  and  down-regulated  genes)  were  enriched  in  this 
function  and  with  more  significant  p-values  at  24  hrs  than  at  1.5  and  4  wks.  Because  the  infil¬ 
trates  only  contribute  transcripts  to  the  tissue,  we  also  analyzed  up-regulated  transcripts  alone. 
These  up-regulated  transcripts  in  AggE  mouse  liver  were  enriched  in  ‘immune  ceU  trafficking’ 
at  24  hrs,  but  not  at  1.5  or  4  wks. 


Figure  4.  Percent  of  metabolites  in  subpathways  elevated  in  AggE  vs.  Ctrl.  Blue  bars  represent  the 
percentage  of  all  metabolites  in  the  subpathway  that  were  elevated  in  plasma  of  AggE  vs.  Ctrl  mice.  Red 
shading  indicatesp  <  0.05  and  light  orange  shading  indicates  p<  0.1  (exact  binomial  test).  Subpathways 
were  included  in  the  figure  only  if  there  was  a  significant  difference  at  these  p-values  in  one  of  the  four  AggE 
vs.  Ctrl  comparisons. 

doi:10.1371/journal.pone.0117092.g004 
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Figure  5.  Higher  levels  of  lipid  and  lower  levels  of  amino  acid  metabolites  in  plasma  of  AggE  vs.  Ctrl 
mice  at  4  wks  after  AggE-withdrawal.  Shown  are  individual  metabolites  that  were  lower  (green)  or  higher 
(red)  in  AggE  vs.  Ctrl  mice  at  4  wks  post-AggE.  Metabolites  are  shown  for  subpathways  that  had  metabolites 
that  were  coordinately  regulated  (shown  in  Fig.  3)  or  that  had  significantly  different  numbers  of  higher  or  lower 
metabolites  in  AggE  mice  (shown  in  Fig.  4).  The  node  labels  indicate  the  superpathway  designation  of  the 
metabolites:  A,  Amino  acid;  C,  Carbohydrate;  L,  lipid;  E,  energy;  P,  Peptide;  and  X,  xenobiotic.  The 
metabolites  are  distributed  using  a  force-directed,  edge-unweighted  algorithm  in  Cytoscape  [89],  with 
biochemical  relationships  shown  as  pink  edges  (KEGG)  and  structural  similarities  shown  as  violet  edges 
(Tanimoto  similarity  >0.7).  Metabolites  for  which  a  KEGG  or  PubChem  Compound  (CID)  Identifier  was 
available  are  shown.  Numbers  in  parentheses  indicate  the  number  of  elevated  (red)  and  suppressed  (green) 
metabolites  in  the  subpathway,  regardless  of  the  availability  of  KEGG  or  CID  identifiers.  Subpathway  labels 
are  shown  next  to  the  groupings  (by  eye)  where  they  were  most  abundant. 

doi:10.1371/journal.pone.0117092.g005 


Upstream  regulator  analysis  also  indicated  inflammatory  changes  of  liver  at  24  hrs,  as  the 
cytokines  ILIB  and  TNF  were  predicted  to  be  activated  with  z-scores  greater  than  2.0  (z-scores 
were  3.2  and  2.2;  andp-values  2.04E-07  and  1.05E-06,  respectively)  [58].  Conversely  at  4  wks 
these  predicted  cytokines  trended  toward  being  inhibited  (ILIB:  z-score  -1.170,  p-value  8.72E- 
04;  TNE:  z-score  -0.303,  p-value  1.77E-05).  Of  all  the  cytokines  with  a  significantly  predicted 
state  of  upstream  activation  for  both  24  hr  and  4  wk  time  points,  16  were  predicted  to  be  in  an 
activated  state  (CCL2,  CSF2,  EDNl,  IFNBl,  ILIA,  ILIB,  IL17C,  IL18,  IL22,  IL23A,  IL27,  IL33, 
LIF,  OSM,  TNF,  and  TNFSF12)  and  2  were  predicted  to  be  in  an  inhibited  state  (IL3, 
TNFSF12).  All  16  with  an  activated  state  were  predicted  for  the  24  hr  samples,  and  the  two  acti¬ 
vators  with  an  inhibited  state  were  predicted  for  the  4  wk  samples. 

Of  all  upstream  activators  at  any  time  point,  the  topmost  one  predicted  for  the  AggE  vs.  Ctrl 
DEGs  was  lipopolysaccharide’  (smallest  p-value)  and  it  was  predicted  to  be  highly  activated  at 
the  24  hr  points  (z-scores  were  4.4  and  3.1;  p-values  5.42E-06  and  4.5E-03,  respectively,  for 
post-5-  and  10-day  stressors),  but  trended  toward  being  inhibited  at  1.5  and  4  wks  (z-scores  -1.9 
and  -1.7;  p-values  2.83E-01  and  9.70E-06,  respectively,  for  5-  and  10-day  stressors). 

At  24  hrs  post-stressor,  despite  no  detectable  difference  in  liver  lymphoplasmacytic  infil¬ 
trates  between  AggE  and  Ctrl  mice,  the  gene  set  enrichment  analysis  indicated  that  the  AggE 
livers  were  more  inflamed  than  the  Ctrl  livers.  Additionally,  well  known  stress  response  genes 
were  among  the  highest  up-regulated  genes  in  AggE  mouse  liver  at  this  acute  time  point  (e.g., 
lipocalin  2,  metallothioneins,  serum  amyloid  proteins).  These  data  suggest  sub-histopatholog- 
ical  liver  damage  was  an  early  effect  of  the  AggE  stress. 

In  contrast  to  the  indicated  inflammation  at  24  hrs  post-stressor,  a  different  picture  emerged 
at  4  wks  post-stressor.  Downstream  function  analysis  in  the  ‘nutritional  disease’  category 
showed  significant  enrichment  in  ‘adult-onset  obesity’.  The  42  enriching  genes  in  ‘nutritional 
disease’  included  26  that  were  enriched  in  ‘weight  gain’.  This  is  consistent  with  the  body  weight 
gain  of  AggE  vs.  Ctrl  mice  that  steadily  increased  up  to  the  10*  day  of  the  stressor  regimen 
(data  not  shown).  In  the  category  ‘hepatic  system  development  and  function’,  significant 
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Table  3.  Topmost  10  canonical  pathway  enrichments  at  4  wks  post-stressor. 


Canonical  Pathway  Name  (genes  enriched)* 

Predicted  downstream  effects 

LPS/IL-1  Mediated  Inhibition  of  RXR  Function 

(ALDH9A1 ,  CYP4A1 1 ,  GSTP1 ,  GSTT1 ,  JUN, 
NR0B2,  SLC27A3,  SREBF1) 

Inhibition  of  cholesterol  and  lipid  metabolism  and 
transport;  activation  of  Phase  1  metabolizing 
enzymes->lipid  and  xenobiotic  metabolism 

PPARa/RXRa  Activation  (ACAA1,  CHUK, 
CYP2C18,  GPD1,  IRS1,  JUN,  NR0B2,  TGFBR1) 

Activation  of  fatty  acid  oxidation,  fatty  acid  uptake, 
lipoprotein  metabolism,  mitochondrial  and 
peroxisomal  beta-oxidation. 

Oxidative  Phosphorylation  &  Mitochondrial 
Dysfunction  (ATP5B,  ATP5D,  COX6A2, 

NDUFA10,  NDUFB5,  SDHA,  SDHB,  UQCRC1) 

Activation 

NRF2-mediated  Oxidative  Stress  Response 

(ACTA2,  EPHX1,  FTL,  GCLM,  GSTP1,  GSTT1, 

JUN) 

Activation  of  phase  1  and  II  metabolizing  enzymes 
and  transport  of  xenobiotics  and  metabolites 

FXR/RXR  Activation  (CYP19A1,  CYP8B1, 

FOXA2,  NR0B2,  PPARG,  SDC1,  SREBF1) 

Inhibition  of  gluconeogenesis,  increased  lipogenesis 
and  lipid  metabolism 

Glucocorticoid  Receptor  Signaling  (CCL1 1 , 
CEBPA,  CHUK,  CSN2,  JUN,  TGFBR1) 

Inconclusive 

Aryl  Hydrocarbon  Receptor  Signaling  (ALDH9A1 , 
CDKN1 B,  GSTP1 ,  GSTT1 ,  JUN,  NR0B2) 

Activation  of  phase  1  and  phase  II  metabolizing 
enzymes 

Superpathway  of  Cholesterol  Biosynthesis 

(DHCR7,  EBP,  FDPS,  HADHA,  NSDHL) 

Activation  of  cholesterol  synthesis 

Type  II  Diabetes  Mellitus  Signaling  (CHUK,  IRS1 , 
PPARG,  SLC27A3,  SLC2A2) 

Reduced  glucose  entry  into  hepatocytes. 

*  Upstream  activators  {p<0.05)  that  additionally  had  a  predicted  activated  activation  state  (z-score  >2) 
totaled  to  12.  These  activators  and  their  target  molecules  totaled  to  109  molecules  and  were  analyzed  for 

overlaps  with  the  shown  canonical  pathways. 

doi:  1 0. 1 371  /journal  .pone.01 1 7092.1003 


enrichment  was  found  in  ‘abnormal  morphology  of  liver’  (liver,  parenchyma,  and  hepatocytes; 
B-H  p-value  =  4.49E-02,  24  genes).  In  the  ‘hepatic  system  development  and  function’  category 
overall,  38  genes  were  enriched.  The  top  six  canonical  pathways  associated  with  these  genes  in¬ 
cluded  ‘PPARa/RXRa  activation’,  ‘aryl  hydrocarbon  receptor  signaling’  and  ‘PPAR  signaling’. 
All  were  predicted  to  be  activated  in  AggE-mice.  Activation  of  the  PPAR  pathways  is  consistent 
with  the  elevated  levels  of  their  ligands  (free  fatty  acids)  in  plasma  [59]  at  4  wks. 

Activation  of  ‘aryl  hydrocarbon  receptor  signaling’  has  been  associated  with  refractoriness 
to  a  second  EPS  challenge,  thought  to  be  part  of  a  “disease  tolerance”  mechanism  [60,61]. 

Thus,  one  possibility  is  that  chronically  activated  ‘aryl  hydrocarbon  receptor  signaling’  could 
have  been  the  consequence  of  liver  exposed  to  gut  bacteria.  Notably,  EPS  (lipopolysaccharide) 
was  the  topmost  predicted  activated  upstream  activator  of  the  AggE-altered  DEGs  at  24  hr 
post-stressor. 

Further  analysis  of  the  4  wk  downstream  enrichments  of  AggE  vs.  Ctrl  DEGs  included  only 
functional  annotations  with  both  a  significant  p-value  for  the  overlap  as  well  as  a  significant  ac¬ 
tivated  state.  This  resulted  in  the  function  annotations:  ‘conversion  of  lipid’,  ‘biosynthesis  of 
nucleoside  triphosphate’,  ‘ascites’,  ‘differentiation  of  granulocytes’,  ‘differentiation  of  myeloid 
cells’,  and  ‘contact  growth  inhibition  of  leukocytes’.  The  immune  cell  functions  involving  dif¬ 
ferentiation  and  growth  inhibition  appear  to  be  consistent  with  the  lymphoplasmacytic  infil¬ 
trates  that  persisted  at  the  chronic  time  points  but  trended  toward  a  lower  grade  at  1.5  and 
4  wks  than  at  24  hrs  (S2  Fig.). 

To  look  further  for  pathways  altered  chronically  by  AggE,  we  performed  upstream  regulator 
analysis  on  the  4-wk  DEGs,  focusing  on  upstream  regulators  with  both  significant  p-values  and 
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z-scores.  Four  were  predicted  to  have  an  inhibited  state  in  AggE  mouse  liver:  RARG  (CD437), 
Map4k4,  HOXC8,  and  D-glucose.  Conversely,  the  activators  with  a  predicted  activated  state  to¬ 
taled  to  12:  SREBFl,  SREBF2,  ESRRA,  ESRRG,  CHUK,  phenobarbital,  TP53.  SCAR,  PPARG, 
mono-(2-ethylhexyl)phthalate,  sulforafan,  and  TO-901317.  Several  of  these  are  plausibly  relat¬ 
ed  through  regulatory  pathways  (e.g.,  mono-(2-ethylhexyl)phthalate,  PPARG,  PPARA,  SCAP, 
SREBFl,  and  SREBF2;  all  found  using  Ingenuity  Pathway  Analysis  (IPA)).  The  12  activated 
upstream  regulators  (excluding  p53)  plus  their  target  molecules  summed  to  109  molecules 
(S3  File).  Analysis  of  their  overlap  with  canonical  pathways  showed  that  the  top  10  canonical 
pathways  (Table  3)  generally  predicted  activation  of  lipoprotein  metabolism,  fatty  acid  uptake 
and  metabolism,  lipogenesis,  and  phase  I  and  II  metabolism  of  lipids  and  xenobiotics;  and  inhi¬ 
bition  of  cholesterol  and  bUe  acid  synthesis  and  transport,  as  well  as  glucose  uptake 
and  gluconeogenesis. 

Consistent  with  our  data,  chronic  psychosocial  stress  has  been  show  to  induce  lipid  dysregu- 
lation,  intrahepatic  accumulation  of  triglycerides,  indicators  of  metabolic  syndrome  [62,63]  as 
weU  as  hepatic  oxidative  stress  and  inflammation  [64] . 

Limitations  in  the  Research  Design 

Our  unbiased  metabolomics  data  showed  high  cj-values  for  pairwise  comparisons  of  metabolite 
levels  between  Ctrl  and  AggE  is  likely  due  to  the  small  sample  size  of  5  or  6  mice  per  group  and 
the  known  highly  variable  stress-response  of  even  inbred  mice.  To  overcome  this  shortcoming, 
we  applied  a  battery  of  statistical  tests,  including  analyzing  metabolite  groupings  constructed 
based  on  KEGG  and  HMDB,  and  the  internal  knowledge  base  of  Metabolon,  Inc. 

Conclusions  and  Hypothesis 

The  repeated  stress  for  5  or  10  days  altered  multiple  organ  systems  acutely,  and  chronic  sys¬ 
temic  changes  were  also  evident.  Acutely  at  24  hrs  post-stressor,  heart  (frank  lesions)  and  liver 
(transcriptomic  changes)  tissues  showed  evidence  of  inflammation,  plasma  protein  levels  indi¬ 
cated  tissue  damage,  and  the  plasma  levels  of  metabolites  derived  from  the  gut  microbiota  were 
altered.  Chronically  at  4  wks  post-stressor,  levels  of  gut-derived  metabolites  were  altered,  hy¬ 
perlipidemia  was  evident,  and  liver  transcriptomics  indicated  activation  of  xenobiotic  and  lipid 
metabolism.  These  data  provide  clues  for  ongoing  efforts  to  understand  and  identify  molecular 
indicators  of  physiologic  responses  to  acute  and  chronic  stressors  that  have  potential  to  predict 
probable  PTSD  and  monitor  its  persistence. 

The  altered  plasma  levels  of  gut-microbiota-derived  metabolites  in  AggE  mice  may  have  re¬ 
sulted  from  the  well-known  effects  of  stress  mediators  on  the  GI  tract  to  induce  changed  motor 
activities,  ulcerations  and  enteritis  [45,46]  and  changed  gut  microbiota  composition  [48,65]. 
Prolonged  stress  has  been  shown  to  produce  low-grade  intestinal  inflammation  [66],  and  in¬ 
flammation  can  increase  translocation  of  commensal  bacteria  across  the  epithelium  that  escape 
from  local  mesenteric  lymph  nodes,  thus  reaching  the  liver  [49].  The  liver’s  ability  to  act  as  a 
firewall  against  systemic  exposure  to  gut  bacteria  could  be  compromised  by  intense  stress,  and 
adrenergic  receptor  activation  is  known  to  constrict  the  hepatic  artery  and  microvasculature, 
reducing  total  hepatic  blood  volume  [67].  Social  disruption  stress  has  been  reported  to  induce 
translocation  of  bacteria  from  the  lumen  of  the  intestines  to  lymphoid  tissue  of  the  body  interi¬ 
or,  and  to  increase  the  percent  of  livers  with  viable  bacteria  [68]  as  well  as  circulating  cytokines 
[46,48].  Such  stress-induced  alterations  could  have  contributed  to  the  altered  liver  gene  expres¬ 
sion  we  observed. 

Stress  mediators  simultaneously  affect  multiple  organ  systems.  Simultaneous  perturbations 
of  closely  cooperating  systems  like  gut  and  liver  could  synergize  to  produce  organ  injury  which 
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leads  to  chronic  changes.  For  example,  liver  disease  resulting  from  bile  duct  ligation  can  cause 
low-grade  intestinal  inflammation  and  histological  damage  to  intestine,  increasing  the  proba¬ 
bility  of  intestinal  microbes  translocating  into  the  vasculature  as  well  as  reducing  the  ability  of 
liver  to  clear  them  [49].  Also,  patients  with  mostly  asymptomatic  liver  steatosis  or  nonalcoholic 
steatohepatitis  (NASH)  had  increased  serum  IgG  and  IgA  responses  against  aerobic  and  anaer¬ 
obic  intestinal  commensal  bacteria.  Additionally,  these  antibody  responses  discriminated  the 
diagnoses  from  healthy  controls  [49] .  Lastly,  dysbiosis  of  intestinal  microbes  has  been  shown 
to  precede  or  trigger  liver  damage  [49,69,70]. 

The  above  findings  collectively  suggest  the  hypothesis  that  repeated  stress  effects  on  multi¬ 
ple  organ  systems  can  synergize  to  result  in  chronically  altered  gut-microbiota-derived  metabo¬ 
lites,  hyperlipidemia,  and  liver  gene  expression  changes.  Such  interactions  of  stress  effects 
among  multiple  organ  systems  might  underlie  reports  of  associations  of  PSTD  with  gastroin¬ 
testinal  complaints,  metabolic  syndrome,  hypertension,  hyperlipidemia,  obesity,  and  cardio¬ 
vascular  disease  [71-75]. 


Methods 

Ethics  statement 

Research  was  conducted  in  compliance  with  the  Animal  Welfare  Act,  and  other  Federal  stat¬ 
utes  and  regulations  relating  to  animals  and  experiments  involving  animals  and  adheres  to 
principles  stated  in  the  Guide  for  the  Care  and  Use  of  Laboratory  Animals  (NRG  2011)  in  facil¬ 
ities  that  are  fully  accredited  by  the  Association  for  the  Assessment  and  Accreditation  of  Labo¬ 
ratory  Animal  Care,  International.  The  protocol  was  approved  by  the  Institutional  Animal 
Care  and  Use  Committee  (lACUC)  of  Walter  Reed  Army  Institute  of  Research,  Silver  Spring, 
MD,  MedStar  Health  Research  Institute,  Washington  DC  and  Animal  Care  and  Use  Review 
Office  (ACURO),  US  Army  Medical  Research  and  Materiel  Command  (USAMRMC). 


Mice 

All  animal  experiments  were  approved  in  writing  by  the  Institutional  Animal  Care  and  Use 
Committee  (lACUC)  of  the  United  States  Army  Medical  Research  and  Materiel  Command 
and  the  Medstar  Research  Institute  (Protocol  Number:  2011-010)  and  were  conducted  in  com¬ 
pliance  with  the  Animal  Welfare  Act,  and  other  Federal  statutes  and  regulations  relating  to  ani¬ 
mals  and  experiments  involving  animals,  adhering  to  principles  stated  in  the  Guide  for  the 
Care  and  Use  of  Laboratory  Animals  (NRC  20 1 1 )  in  facilities  fully  accredited  by  the  Associa¬ 
tion  for  the  Assessment  and  Accreditation  of  Laboratory  Animal  Care,  International.  All  mice 
were  purchased  from  Jackson  Laboratory,  Bar  Harbor,  ME,  USA.  AU  mice  had  free  access  to 
food  and  water  and  were  kept  in  a  temperature-controlled  room  (21  -i-  2°C)  on  reverse  12/12  h 
light/dark  cycle  (lights  off  at  06:00  AM). 

Aggressor  mice.  The  SJL  albino  male  mice  (5  to  6  weeks  old  and  weighing  30-35  g  when 
purchased),  which  are  characterized  by  extreme  aggression  in  males  [76],  were  housed  individ¬ 
ually  in  polycarbonate  cages  (48  X  27  X  20  cm)  for  one  month  prior  to  the  experiment  to  in¬ 
duce  aggressiveness  due  to  isolation.  They  were  then  trained  to  attack  intruders  as  described 
[11]. 

Subject  mice.  Male  C57BL/6J  mice  (5  to  6  weeks  old  weighing  20-25  g)  were  singly  housed, 
in  a  different  room  from  the  aggressor  mice  under  the  same  environmental  conditions,  for  one 
week  prior  to  and  during  the  stress  regimens,  and  afterward  during  the  1.5  or  4  wk  stress- 
withdrawal  periods. 

Control  mice.  Control  mice  were  housed  identically  to  subject  mice. 
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Aggressor  exposure 

Aggressor  exposure  (AggE)  sessions  followed  a  modified  “cage-within-cage  resident-intruder” 
protocol.  For  aggressor  exposures,  subject  mice  were  placed  in  a  wire  mesh  cage  (17.5  X  14  X  7.5 
cm)  that  was  put  inside  the  aggressor 's  large  plastic  home  cage  (48  X  27  X  20  cm)  for  6  hrs  as 
described  [11].  During  the  6  hrs,  mice  were  randomly  placed  in  physical  contact  with  the  aggres¬ 
sor  mouse  3  times,  each  for  1  min.  or  10  strikes  (whichever  came  first).  The  6-hr  sessions  were 
repeated  daily  for  either  5  or  10  consecutive  days.  In  a  separate  room,  control  mice  were  sham- 
AggE  using  the  same  cage-within-cage  environment  and  food  deprivation  on  the  same  6-hr 
daily  regimen,  but  with  fresh  bedding  and  without  any  aggressor  present.  Aggressor  mice  were 
provided  food  and  water  ad  libitum  during  the  6-hr  “cage-within-cage”  sessions.  After  each  6-hr 
session,  AggE  and  sham- AggE  (Ctrl)  mice  were  single-housed  in  their  home  cages  with  ad 
libitum  food  and  water  until  the  next  session.  Twenty-four  hours  after  the  5  and  10  days  of 
daily  6-hr  sessions,  half  of  each  population  was  randomly  chosen  and  Ctrl  and  AggE  were 
euthanized  alternately.  The  other  half  of  the  population  continued  to  be  single-housed  in 
their  home  cages  until  they  were  euthanized  at  1.5  weeks  (5-day  AggE)  or  4  weeks  (10-day 
AggE). 

One  hour  prior  to  the  terminal  bleed,  mice  were  subjected  to  a  partition  test  to  evaluate  be¬ 
haviors  in  response  to  the  situational  reminder  of  the  aggressor’s  presence.  Each  subject  mouse 
was  placed  into  an  aggressor’s  home-cage  on  the  other  side  of  a  partition  that  blocks  physical 
interaction  between  the  mice  but  permits  passage  of  sensory  cues.  Subject  mouse  were  then 
video  recorded  for  five  minutes  for  behavioral  evaluations  [11]. 


Tissue  collection  and  processing 

Mice  were  euthanized  by  cervical  dislocation.  Ctrl  and  AggE  mice  were  alternately  euthanized 
and  processed  for  tissue  collection  immediately  after  euthanasia.  Plasma,  liver  and  duodenum 
were  frozen  immediately  until  analysis.  Heart,  stomach,  liver,  gallbladder,  and  kidney  were  pro¬ 
cessed  for  histopathology. 

Blood  was  collected  from  heart  with  a  20  gauge  needle  first  rinsed  with  sodium  citrate,  and 
was  transferred  to  tubes  containing  50  pi  of  0.015  M  buffered  sodium  citrate.  Tubes  were 
mixed  gently  three  times  and  then  centrifuged.  Plasma  was  transferred  to  a  new  tube  and  cen¬ 
trifuged  again.  The  resultant  clean  and  clear  plasma  was  collected,  immediately  placed  on  dry 
ice,  kept  frozen,  and  shipped  frozen  to  Metabolon,  Inc.  (Durham,  NC),  where  the  frozen  plas¬ 
ma  was  immediately  stored  at  -80°C. 

For  liver  gene  expression  analysis,  liver  was  weighed  before  thawing  and  Trizol™  was 
added  prior  to  homogenization  and  RNA  isolation  according  to  the  manufacturer’s  recom¬ 
mendations  (Invitrogen,  Grand  Island,  NY).  The  RNA  was  cleaned  using  the  RNAeasy  mini 
column  (Qiagen).  Extracted  RNA  from  all  samples  was  stored  -80°C.  RNA  quality  and  quantity 
were  assayed  using  a  Nanodrop  ND-1000  spectrometer  (Thermo  Fisher,  Wilmington,  DE). 
RNA  from  each  sample  was  assessed  for  purity  by  A260/280  ratios.  To  determine  the  integrity 
of  rRNA,  1  pL  of  eluate  from  each  sample  was  analyzed  using  the  Agilent  2100  Bioanalyzer  or 
Agilent  2200  Tapestation  (Agilent  Technologies,  Inc.,  Santa  Clara,  CA)  according  to  the 
manufacturer’s  instructions. 


Plasma  metabolomics 

Unbiased  metabolomic  profiles  of  plasma  samples  were  obtained  using  ultrahigh  performance 
liquid  chromatography/tandem  mass  spectrometry  (UHPLC/MS/MS)  and  gas  chromatogra¬ 
phy/  mass  spectrometry  (GC/MS  [77,78].  MS/MS  data  were  searched  against  a  library  of 
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purified  standards  or  recurrent  unknown  entities  [77],  and  the  resulting  named  metabolites 
were  used  for  our  analysis. 

Samples  were  extracted  to  remove  protein  and  dislodge  small  molecules  bound  to  protein  or 
physically  trapped  in  the  precipitated  protein  matrix.  Extracts  were  split  for  analysis  on  the  gas 
chromatography  (GC)  and  liquid  chromatography  (LC)  platforms. 

Sample  extracts  for  LC/MS  (LC/MS,  LC/MS/MS)  were  split  into  two  aliquots,  dried,  and 
then  reconstituted  in  acidic  or  basic  LC-compatible  solvents,  each  of  which  contained  1 1  or 
more  injection  standards  at  fixed  concentrations.  One  aliquot  was  analyzed  using  acidic  posi¬ 
tive  ion  optimized  conditions  and  the  other  using  basic  negative  ion  optimized  conditions  in 
two  independent  injections  using  separate  dedicated  columns.  Extracts  reconstituted  in  acidic 
conditions  were  gradient  eluted  using  water  and  methanol  (both  containing  0.1%  formic  acid) 
while  the  basic  extracts,  which  also  used  water/methanol,  contained  6.5  mM  ammonium  bicar¬ 
bonate.  The  MS  analysis  alternated  between  MS  and  data-dependent  MS^  scans  using  dynamic 
exclusion.  The  system  was  based  on  a  Waters  ACQUITY  UPLC  and  a  Thermo-Einnigan  LTQ 
mass  spectrometer,  which  consisted  of  an  electrospray  ionization  (ESI)  source  and  linear  ion- 
trap  (LIT)  mass  analyzer. 

Samples  destined  for  GC/MS  analysis  were  re-dried  under  vacuum  desiccation  for  a  mini¬ 
mum  of  24  hs  prior  to  being  derivatized  under  nitrogen  using  bistrimethyl-sUyl-triflouroaceta- 
mide  (BSTFA).  The  GC  column  was  5%  phenyl  and  the  temperature  ramp  was  from  40°  to 
300°C  during  a  16  minute  period.  Samples  were  analyzed  on  a  Thermo-Einnigan  Trace  DSQ 
fast-scanning  single-quadrupole  mass  spectrometer  using  electron  impact  ionization.  The  in¬ 
strument  was  tuned  and  calibrated  for  mass  resolution  and  mass  accuracy  on  a  daily  basis.  The 
information  output  from  the  raw  data  files  was  automatically  extracted. 

For  Accurate  Mass  Determination  and  MS/MS  fragmentation  (LC/MS),  (LC/MS/MS),  the 
LC/MS  was  based  on  a  Waters  ACQUITY  UPLC  and  a  Thermo-Finnigan  LTQ-FT  mass  spec¬ 
trometer,  which  had  a  linear  ion-trap  (LIT)  front  end  and  a  Fourier  transform  ion  cyclotron 
resonance  (FT-ICR)  mass  spectrometer  backend.  For  ions  with  counts  greater  than  two  mil¬ 
lion,  an  accurate  mass  measurement  could  be  made  on  the  parent  ion  as  well  as  fragments.  The 
typical  mass  error  was  less  than  5  ppm.  Ions  with  less  than  two  million  counts  require  a  greater 
amount  of  effort  to  characterize.  Fragmentation  spectra  (MS/MS)  were  typically  generated  in 
data  dependent  manner.  Metabolon’s  proprietary  visualization  and  interpretation  software  was 
used  to  confirm  the  consistency  of  peak  identification  among  samples. 

Compounds  were  identified  by  comparison  to  purified  standards  or  recurrent  unknown 
entities.  The  purified  standards  consisted  of  more  than  1800  commercially  available  com¬ 
pounds.  Match  to  the  specific  compound  or  an  isobaric  entity  was  indicated  by  the  combina¬ 
tion  of  chromatographic  properties  and  mass  spectra.  Additional  entities  were  identified 
as  recurrent  entities  (both  chromatographic  and  mass  spectral)  and  currently  remain 
unidentified. 

Superpathways  and  subpathways  were  constructed  from  information  available  in  public  da¬ 
tabases  (e.g.,  KEGG  and  HMDB  and  textbooks)  and  Metabolon’s  internal  knowledge  base. 
They  generally  reflect  structural  and  substrate-product  relationships  and/or  the  biology  of  the 
gut  microbiota.  These  superpathways  and  subpathways  were  constructed  prior  to  the 
statistical  analyses. 

Corticosterone  ELISA.  Blood  collected  in  tubes  containing  50  pi  0.105  M  buffered  sodium 
citrate  was  centrifuged  to  separate  plasma  that  was  stored  at  -80°C  until  the  assay.  The  ELISA 
was  based  on  a  competitive  enzyme  immunoassay  using  a  combination  of  specific  antibodies 
to  corticosterone  and  HRP-labeled  corticosterone  (Kamiya  Biomedical  Company,  Seattle, 

WA). 
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Plasma  Proteins 

Blood  samples  were  collected  in  plastic  tubes  containing  buffered  sodium  citrate  (3.2%)  and 
were  centrifuged  to  separate  plasma.  Seventy  pi  of  samples  were  shipped  to  Rules-Based  Medi¬ 
cine,  Inc.  (Austin,  TX).  All  samples  were  stored  at  -80°C  until  tested  and  were  analyzed  for  Ro¬ 
dent  Multi- Analyte  Profile  (MAP)  v.2.0  Antigens  (Myriad- RBM)  using  Luminex  technology. 

Liver  microarray  gene  expression  analysis 

mRNA  was  extracted  from  liver  samples  of  4  to  6  mice  from  each  of  the  eight  different  time/ 
treatment  groups  and  were  processed  for  dual  dye  microarrays  (Whole  Mouse  Genome  Micro¬ 
array  Kit,  Agilent  Technologies,  Inc.,  CA)  following  the  vendor’s  protocol.  Purified  RNA  was 
labeled  with  Cy-5  and  the  reference  RNA  (Agilent  Technologies,  Inc.)  with  Cy-3.  The  samples 
were  simultaneously  hybridized  to  Agilent  4  x  44k  slides  (platform  number:  14868)  for  17  h  at 
55°C.  After  hybridization,  slides  were  processed  in  a  series  of  washes  and  scanned  using  an  Agi¬ 
lent  DNA  microarray  scanner.  Features  were  extracted  using  the  default  setting  of  the  Feature 
Extraction  software  v.10.7  (Agilent  Technologies,  Inc.).  The  resulting  expression  data  files  were 
analyzed  using  custom  R  scripts.  Arrays  were  background  corrected,  individually  LOWESS- 
normalized  [79],  and  quantile-normalized  across  all  arrays.  Probes  were  then  filtered  to  include 
only  those  associated  with  a  REFSEQ  gene  symbol  (~26.5k  out  of  45k).  The  remaining  probes 
were  compared  between  groups  using  a  moderated  t-statistic  with  Benjamini-Hochberg  false 
discovery  rate  (FDR)  correction.  The  microarray  data  was  submitted  to  the  Gene  Expression 
Omnibus  (GEO).  This  can  be  searched  using  the  Platform  ID:  GPL4134,  Series:  GSE58275. 

Quantitative  Real-time  PCR 

Microarray  results  were  validated  using  qPCR  primer  assays  (SABiosciences,  Frederick,  USA) 
using  Applied  Biosystems  7900  systems.  The  RNA  of  each  mouse  was  converted  into  cDNA 
using  the  RT^  First  Strand  Kit  (SABiosciences,  Frederick,  USA).  This  cDNA  was  then  added  to 
the  RT^  SYBR  Green  qPCR  Master  Mix  (SABiosciences,  Frederick,  USA).  All  steps  were  done 
according  to  the  manufacturer’s  guidelines.  The  data  were  analyzed  using  AAct  method  [80] 
using  18S  ribosomal  RNA,  fiactin  and  Ubiquitin  C  as  housekeeping  genes.  Validation  of  micro¬ 
array  data  for  selected  genes  is  shown  in  S3  Fig. 

Histopathology 

Organ  samples  were  collected  from  euthanized  control  and  AggE  mice  and  perfused  in  ice- 
cold  4%  paraformaldehyde.  Tissues  were  then  stained,  sectioned  to  5  pm  thickness,  and 
mounted.  A  board-certified  veterinary  pathologist,  blinded  from  the  animal  condition,  identi¬ 
fied  and  analyzed  the  samples  using  brightfield  optics.  Histologic  findings  were  scored  as  nor¬ 
mal,  minimal,  mdd,  moderate,  and  marked,  as  well  as  focal,  multifocal  and  diffuse. 

Statistical  analysis 

Levels  of  metabolites  or  proteins  in  plasma  between  mouse  groups  were  compared  using  pair¬ 
wise  methods  including  Welch’s  t-test  and  ANOVA.  Random  forests  analysis  was  used  to  iden¬ 
tify  metabolites  that  classify  Ctrl  and  AggE  mouse  groups.  The  groups  receiving  the  5 -day  and 
10-day  stressors  were  analyzed  individually  and  combined.  We  also  used  single-node  decision 
trees  (Gini  index  and  Information  Gain  (KuUback-Leibler  divergence)).  Chi-square  test,  Kol¬ 
mogorov- Smirnov  test,  and  an  in-house  procedure  that  examines  samples  with  near-extremum 
intensities,  which  is  a  simple  classifier  that  uses  a  leave-one-out  nearest-centroid  procedure  for 
classification.  The  five  top  scoring  metabolites  from  the  methods  were  selected  as  candidate 


PLOS  ONE  I  DOI:10.1371/]oumal.pone.01 17092  January  28,  2015 


18/24 


PLOS 


ONE 


Metabolomics  and  Liver  Transcriptomics  in  PTSD  Mouse  Model 


discriminators  and  were  correlated  with  each  other.  Those  that  correlated  with  a  Pearson’s  cor¬ 
relation  coefficient  >  0.7  were  placed  within  one  category,  and  only  the  metabolite  within  each 
category  with  the  largest  maximum  intensity  was  evaluated  further.  These  metabolites  were 
correlated  with  all  other  metabolites  to  identify  those  potentially  belonging  to  related  metabo¬ 
lite  pathways  (Pearson’s  correlation  coefficient  >  0.7,  p-value  of  0.024  for  any  comparison).  All 
metabolites  that  were  exact  separators  by  Gini  index  and  the  remaining  metabolites,  plus  those 
that  correlated  with  them,  were  used  in  exhaustive  searches  with  distance-dependent  K-Nearest 
Neighbor  (DD-KNN)  and  rule-based  classifiers  to  identify  sets  of  metabolites  that  best  discrim¬ 
inate  Ctrl  from  AggE  mice. 

Comparisons  of  up  vs.  down  metabolites  in  Ctrl  vs.  AggE  at  each  time  point  were  done 
using  a  two  tailed  binomial  exact  test  against  the  null  hypothesis  of  a  50%  up  and  50%  down,  as 
implemented  in  the  R  binom.test  function  [81-83].  Comparisons  of  differences  in  the  propor¬ 
tion  of  up  vs.  down  metabolites  between  time  points  were  done  using  a  two-tailed  difference  of 
proportions  test  against  the  null  hypothesis  of  equal  proportions,  as  implemented  in  the  R 
prop.test  function  [83-87]. 

A  Manhattan  distance  method  was  used  to  identify  subpathways  having  metabolites  that 
were  altered  similarly  among  the  eight  simple  comparisons  of  AggE  vs.  Ctrl  and  24  hrs  vs.  1.5 
or  4  wks.  Each  metabolite  of  a  subpathway  was  assigned  either  a  0  or  a  1  for  having  its  mean 
plasma  level  higher  or  lower  in  each  comparison.  Here,  the  actual  direction  of  regulation  is  not 
important,  only  the  pattern  relative  to  other  metabolites.  To  get  the  average  distance  of  metab¬ 
olites  in  a  subpathway,  each  metabolite  of  a  subpathway  was  compared  pairwise  with  every 
other  metabolite  of  the  pathway.  Two  metabolites  with  the  same  Ts  and  O’s  in  all  eight  compar¬ 
isons  (same  pattern  of  regulation  across  experimental  conditions)  have  a  distance  of  zero.  Two 
metabolites  with  different  Ts  and  O’s  in  all  eight  comparisons  have  a  distance  of  8.  The  total  up 
vs.  down  Manhattan  distance  across  all  pairs  of  metabolites  within  a  subpathway  was  tested 
against  the  null  distribution  of  analogous  distances  created  by  taking  10,000  random  samplings 
from  aU  the  metabolites  of  the  number  in  the  subpathway  of  interest.  An  approximate  (slightly 
conservative)  p-value  was  then  calculated  as  (b+l)/(m+l),  where  b  is  the  number  of  sample  dis¬ 
tances  greater  than  the  actual  subpathway  distance  and  m  is  the  number  of  random  samplings 
(10,000),  as  described  in  [88]. 

Gene  Enrichment  Analysis.  Data  were  analyzed  using  causal  analysis  approaches  in  Inge¬ 
nuity  Pathway  Analysis  (IPA,  Ingenuity  Systems,  Inc.,  Redwood,  CA)  to  infer  upstream  biolog¬ 
ical  causes  and  probable  downstream  biological  effects,  using  enrichment  score  (Eisher’s  Exact 
Testp-value)  to  measure  overlap  between  observed  and  predicted  gene  sets,  and  z-score  to 
match  observed  and  predicted  up/down  regulation  patterns  [58]. 


Supporting  Information 

51  Fig.  Cardiac  histopathology  of  AggE  mice.  Histopathology  was  scored  as  normal,  mini¬ 
mal,  mUd,  moderate,  as  well  as  focal,  multifocal  and  diffuse.  Number  of  mice  evaluated: 

5  Day-24  hrs:  Ctrl  =  7;  AggE  =  8.  5  Day- 1.5  wks:  Ctrl  =  5;  AggE  =  5.  10  Day-24  hrs:  Ctrl  =11; 
AggE  =11.  10  Day-4  wks:  Ctrl  =  5;  AggE  =  5. 

(TIFE) 

52  Fig.  Histology  of  liver,  gallbladder,  stomach  and  kidney.  Histology  was  scored  as  normal, 
minimal,  mUd,  moderate,  as  well  as  focal,  multifocal  and  diffuse.  Number  of  mice  evaluated: 

5  Day-24  hrs:  Ctrl  =  7;  AggE  =  8.  5  Day- 1.5  wks:  Ctrl  =  5;  AggE  =  5.  10  Day-24  hrs:  Ctrl  =11; 
AggE  =11.  10  Day-4  wks:  Ctrl  =  5;  AggE  =  5. 

(TIFF) 
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S3  Fig.  Validation  of  gene  expression  changes  observed  in  microarray  experiments. 

(TIFF) 

51  File.  Unbiased  plasma  metabolomics  analysis. 

(XLSX) 

52  File.  Analysis  of  plasma  proteins. 

(XLSX) 

53  File.  Upstream  regulators  plus  their  target  molecules  at  4  wks. 

(XLSX) 

51  Table.  Summary  of  altered  metabolites. 

(TIFF) 

52  Table.  Liver  gene  enrichment  suggests  inflammation  at  24  hrs  post-5-day  AggE  was  re¬ 
solved  at  24  hrs  post- 10-day  AggE. 

(TIFF) 
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